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Abstract 

We present results for the kaon i?-parameter Bk from a sample of 200 configurations 
using the Wilson action and 460 configurations using the SW-Clover action, on a 
18 3 x 64 lattice at (3 = 6.0. We compare results obtained by renormalizing the relevant 
operator with different "boosted" values of the strong coupling constant a s . In the case 
of the SW-Clover action, we also use the operator renormalized non-perturbatively. In 
the Wilson case, we observe a strong dependence of Bk on the prescription adopted for 
a s , contrary to the results of the Clover case which are almost unaffected by the choice 
of the coupling. We also find that the matrix element of the operator renormalized 
non-perturbatively has a better chiral behaviour. This gives us our best estimate of 
the renormalization group invariant i?-parameter, Bk = 0.86 ± 0.15. 
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1 Introduction 



The kaon £?-parameter Bk is a quantity of great phenomenological interest, being 
related to the e-parameter which measures CP-violation in the K°-K° system. With 
a large top quark mass Q , an accurate prediction of Bk enables us to limit the range 
of values of the CP-violation phase 6; see for example ref. Bk is obtained from the 
calculation of the matrix element of the relevant four-fermion operator (K \O AS=2 \K°) , 
where 0^ s=2 = sj^dsj^d; s and d stand for strange and down quarks and 7^ = 
^7^(1 — 75). A non-perturbative estimate of the matrix element may be obtained in 
the framework of lattice QCD ||, |J]. Several lattice calculations of Bk have been 
performed in recent years, both with staggered |5|, |6| and Wilson [j7fl- |lC[] fermions. 

With the Wilson or SW-Clover quark lattice actions, due to the presence of the 
chiral symmetry breaking term, AS=2 mixes with operators of different chirality [11 



12|. For this reason, it is possible to define a renormalized operator with definite 
chiral properties only in the continuum limit, i.e. when a — > 0. At finite a, one 
can improve the chiral behaviour of the matrix element of AS = 2 by subtracting a 
suitable set of dimension-6 operators. The mixing coefficients have been computed 
so far only in one-loop perturbation theory [12] — [ 15] . In this way, the systematic 
error in the value of the matrix element determined on the lattice is of 0(a 2 ) and, 
due to the finitess of the lattice spacing, of 0(a). Following the Symanzik proposal, 
one can reduce the discretization errors from 0(a) to 0(a s a) by using the tree-level 



improved" SW-Clover quark action jig, 17]. Using this action, the improvement has 



been shown to be effective for two-fermion operators, at values of (5 currently used in 
numerical simulations [18]-p0]. It remains true, however, that ignorance of higher- 
order perturbative corrections to the mixing coefficients, which are expected to be 
large [21], can distort the chiral behaviour of the operator and hence induce a large 
systematic error in the determination of Bk- For this reason the kaon matrix element 
of 0^ s=2 does not vanish in the chiral limit 22]. This renders the evaluation of Bk 
more problematic with Wilson or SW-Clover fermions than with staggered fermions. 

In this paper, we present an extended study of (K \O AS=2 \K°) , using different sets 
of quenched gauge field configurations, at f3 = 6.0. We address the two systematic 
effects which afflict its chiral behaviour: corrections of 0(a 2 s ) or higher in the mixing 
coefficients and contributions of 0(a) due to the finiteness of the lattice spacing a. In 
order to estimate both sources of error, we have performed a high statistics analysis 
of these effects. Discretization errors have been studied by comparing the results for 
Bk obtained with the Wilson and SW-Clover actions. The systematic error due to the 
use of perturbation theory in the calculation of the mixing coefficients has also been 
analysed by studying the effects of a different choice of the expansion parameter, as 
proposed in ref. [21]. We also used, in the Clover case, the coefficients computed in 



[23] with the Non-Perturbative Method (NPM) of ref. |24l]. The quality of our results 
is somewhat limited by the use of the "thinning" procedure, imposed by memory 
limitations of the APE 6-Gflops machine pl| ]. Our tentative conclusion is that the 
most appreciable improvement of the chiral behaviour of the matrix element comes 
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from the implementation of the non-perturbative calculation of the mixing coefficients, 
whereas Clover improvement and Boosted Perturbation Theory (BPT) seem to be less 
influential. Preliminary results of this analysis have been presented in |26| . 

Section ^ introduces the notation necessary for the renormalization of the lattice 
operator; in section || we spell out all the details of the numerical simulation; the main 
results of our work are contained in section |I]; in particular the various fitting procedures 
are explained and the effects of Clover improvement, boosting and non-perturbative 
subtraction are analysed; in section [5] we extract the renormalization group invariant 
value of Bk', section ^ contains our conclusions. 

2 The renormalized operator 

The weak effective Hamiltonian relevant to K°~K° mixing is 

n^f 2 = c{M w /^)d^= 2 ^) (i) 

where O as=2 (/j) is the renormalized operator; C{My/ j [i) is the corresponding Wilson 
coefficient in the Operator Product Expansion (OPE) and // is the renormalization 
scale. Our aim is to measure the matrix element (K \O AS=2 (fi)\K°), parametrized by 
the Bk parameter as follows: 

(K»\6* S =\»)\K») = \fW K BK{») , (2) 

where 

(K \d AS = 2 (»)\K°) - ^^ | S - 7A1 d|0}(0| S - 7 ^|^ ) = \f K m\ (3) 

is the matrix element in the vacuum saturation approximation. 

Since the lattice regularization with Wilson fermions explicitly breaks chiral symme- 
try, the lattice AS=2 mixes under renormalization with other "effervescent" operators 
with the same flavour numbers but with different chirality [pjj| -|15]: 



6 AS=2 (n) = Z + (na, a s ) [6 AS = 2 (a) + Z SP 6 SP (a) + Z VA 6 VA (a) + Z SPT d S PT («)] , 

( 4 ) 

where Zgpy a,spt = Z spy a,spt(o> s )] 0(a) denotes the bare lattice operators 

1 

Osp(a) = -j^[sdsd- sj 5 dsj 5 d}] (5) 
iV 2 + 7V-l_ 

Ova (a) = ^ [sj^ds^d - sj^ds-f^d] ; (6) 

N - 1 

Ospt (a) = N [sdsd + sj 5 dsj 5 d + sa^dsa^d] . (7) 

and N is the number of colours. In the chiral limit, the basis of operators appearing 
on the r.h.s. of eq. (0) remains unchanged beyond one-loop, due to CPS symmetry 
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[27|. Note that, since there are no AS = 2 operators of lower dimension, the mixing 
takes place between dimension-6 operators only. This implies that the mixing coeffi- 
cients are at most logarithmically divergent and can be computed on the lattice using 



Perturbation Theory (PT) uM. This calculation has been performed to one-loop for 



both the Wilson and Clover actions [12, 13, 15 1; the results are expressed in terms of 
two constants F + and F* defined as 

Z+ = 1 + ~r-F + , Zsp = Zva = Zspt = -r-F* ' . (8) 

Notice that the equality between the off-diagonal mixing coefficients holds only at 
one-loop. Their values are given in table |l[ In the Wilson case, F + is the constant 
necessary to relate the bare lattice operator to the continuum operator renormalized 
at fi = 1/a in the MS Dimensional Reduction scheme (DRED) [T^] . The perturbative 
result for F + in the Clover case requires some explanation. In ref. |]l^], F+ = —10.9 
is the constant necessary to relate the bare lattice operator to the continuum operator 
renormalized at [i = 1/a in DRED. The value of F+ reported in the table, instead, refers 
to the overall lattice renormalization constant necessary to normalize the operator in 



the Regularization Independent (RI) scheme, defined in ref. H and used in ref. [23 
in the Landau gauge (F + — > F + — 8 In 2 + 5/3 + L(pt)). This is the same scheme used 
in the non-perturbative case, to which we want to compare the perturbative results. 
In the RI scheme, there is a dependence of F + on the momentum p 2 = fi 2 of the 
external states, which is reflected in the term L(fi) = 2/(167r 2 ) ln(fi 2 a 2 ), appearing in 
the table. This term is numerically small in the range of fi 2 a 2 used in our study: it 
gives a contribution of order 0.01 to Z+, if one uses the bare coupling, and of order 
0.02 if one uses some boosted coupling, see below; for this reason, it will be ignored in 
the following. 

Although the renormalized operator 0^ s=2 (fi) is constructed so as to have the 
correct chiral properties in the continuum limit, the perturbative estimate of the Z's 
is bound to introduce a systematic distorsion arising from the one-loop truncation. 
BPT should improve the behaviour of the perturbation series, thus providing better 
estimates of the Z's |pl |. On the other hand, the NPM for the calculation of the 
Z's, as proposed in ref. | 24}| , should in principle account for all perturbative and non- 
perturbative contributions. In order to compare the perturbative mixing coefficients 
to the non-perturbative ones, in table [1] we also give the effective F's obtained non 



perturbatively in the Clover case [23|. Since in the non-perturbative case there is some 
dependence on the renormalization scale [J? a 2 , see also below, we give the results at 
two values of p?a 2 . The F's are defined from Z$py 'A,spt = Qt-s /(^)Fsp,va,spt, with 
the Z's taken from the non-perturbative calculation. We denote by = 1.68 a s the 
coupling constant in the Boosted Scheme (BS) 1, to be introduced in section |], and by 
a s = 3/(271-/3) the bare lattice coupling constant. 

Were the coefficients of the mixing with the dimension 6 operators to be known with 
infinite accuracy, a further source of systematic error would still be present, due to the 



finiteness of the lattice spacing. The SW-Clover action [16] removes O(a) effects from 
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Method 


F + 


^5P 


F V A 


Fspt 


W-PT 


-47.6 


9.6 


9.6 


9.6 


SW-PT 


-14.8 + L(p) 


19.4 


19.4 


19.4 


NPM-fi 2 a 2 = 0.96 


-15(3) 


13(7) 


28(2) 


23(4) 


NPM-/iV = 2.47 


-14(2) 


21(6) 


31(2) 


22(3) 



Table 1: Renormalization constants for the Wilson (W) and SW-Clover (SW) actions. The 
first two rows are the results from one-loop perturbation theory, denoted by F + and F* in 
the text. The last two lines are the results for the effective constants F +j sp,va,spt, obtained 
from the corresponding Z +t spyA,sPT, calculated with the NPM in the Clover case [23], at 
two different values of the renormalization scale. 



all matrix elements, leaving us with 0(a s a) and 0(a 2 ) systematic corrections [lTj. It 
must be noted that in the Clover case all the results refer to the "improved-improved" 



operators introduced in refs. [15 



3 Computational details 

The calculations have been done in the quenched approximation, generating 460 gauge 
configurations on a 18 3 x 64 lattice at [5 = 6.0 with a standard Metropolis algorithm, 
using the 6-Gflops version of APE. The light fermion propagators have been com- 
puted on the full set of configurations using the SW-Clover action, and on a subset 
of 200 configurations using the Wilson action. The reason for this difference is the 
following. Results obtained with the two actions on the same ensemble of 200 con- 
figurations demonstrated that the Wilson action is characterised by smaller statistical 
fluctuations. In order to have comparable errors, we then increased the statistics to 
460 configurations in the Clover case. The values of the hopping parameter used were 
K = 0.1550,0.1540,0.1530 for the Wilson case and K = 0.1440,0.1432,0.1425 for the 
Clover case. All two- and three-point correlation functions have been calculated at 
degenerate quark masses. The statistical error has been estimated by the jacknife 
method, by decimating 20 (46) configurations at a time, for the ensembles of 200 (460) 
configurations. 

The 0^ s=2 operator was placed at the origin. In what follows, 0^ s=2 denotes 
the renormalized operator defined in eq. @. One kaon source was inserted at fixed 
time t y = 12 or 16, whereas the second kaon source was allowed to vary over the full 
time-range t x = 0, 63; each kaon source carries a spatial momentum denoted by p or q. 
In lattice units of 2ir/L (L is the number of sites in each spatial direction) the values 
considered were po = (0, 0,0), p± = (±1, 0, 0),pn = (1,1,0) and the corresponding 
ones equivalent under the cubic group. The same values have also been attributed to 
q. The combinations of momenta examined were (po,qo); (po,Q±)] (Po,Qn)', (p±,Q±)- 
The last case includes three possibilities of relative momentum orientation, namely 
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parallel, antiparallel and perpendicular. Cases which are equivalent under the cubic 
group have been averaged. 

The meson two-point functions were fitted in the time interval t x = 10-22, which 
guarantees a good isolation of the lightest pseudoscalar state. By extrapolating (linear- 
ly in the quark mass) to the vanishing pseudoscalar mass, we computed the critical 
hopping parameter K cr = 0.15697(3) for the Wilson action and K cr = 0.14546(2) for 
the SW-Clover action. The lattice spacing determined from the p mass is a -1 = 2.19(4) 
GeV (W) and a -1 = 2.06(4) GeV (SW). Using the kaon mass and interpolating the 
lattice pseudoscalar masses we find the hopping parameter value at the strange quark 
mass, K s = 0.15468(9) (W) and K s = 0.14366(7) (SW)._ 

In order to extract Bk from the matrix element (K° (p)\0 =2 ([i)\K° (q)) , using 
the method exposed in section ||, we calculated the following two- and three-point 
correlation functions 

G 5 (t x ;p) = ^(s(x) 75 d(x)d(0) 7 5 S (0)) e -^, 

X 

G A (t x ) = ^(s(xho75#c)J(0)75s(0)) , (9) 

X 

G 6 (t x ,t y ;p,q) = ^(d(y) 75 Ky)O A5=2 (0)d>) 75 ^x))e-™ e ^, 



and the ratios 



^ ^ G 6 (t x ,t y -p,q) (K\p)\6^= 2 \K»{q)) ^ 
G 5 (t x ;p)G 5 (t y ;q) Z 5 



_ 8 G\{t x ) ^ 8 /| ^2 
SGj(t x ;0) 3Z\Z 5 J 



where Z§ = (0\s(0)^^d(0)\K }; fx is the kaon decay constant and Za is the axial 
current renormalization constant. 

The asymptotic behaviour indicated by the arrows in the above equations is valid 
in the limit t y <C <C t x . In the absence of final state interaction [29] and finite 



volume effects, the ratio R% with time ordering <C (t y ,t x ) would isolate the matrix 
element (K°(p)K (-q)\6 AS=2 \0). On our lattice with periodic boundary conditions 
both signals appear as two distinct plateaux. In fig. |l|, we show two typical cases 
among the best (fig. [|a) and the worst (fig. |l|b) plateaux for R3. Two plateaux are 
clearly visible; the rightmost one corresponds to the case of interest (^' |O A5=2 |i(' ), 
whereas the leftmost one corresponds to (i^°K°|d A5=2 |0). 

Since one can average those cases characterised by spatial momenta equivalent 
under the cubic symmetry, we are faced with two options: either average the correlation 
functions before taking their ratios or average the ratios. Having found no significant 
difference between the results obtained from the two methods, we have opted for the 
first one. For example, in fig. Rib, we show the ratio of the average of three correlations. 
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Figure 1: R3 for the SW-Clover action at K = 0.1440 (a) p — q — po; (b) average of three 
cases: p '= q '= p + = (1,0, 0); p = q= (0, 1, 0) ; p = q = (0, 0, 1). The renormalized operators 
have been obtained in perturbation theory, by using ctp as the boosted coupling, cf. eq. (17). 



We note that the results obtained from the data with t y = 16 are fully compatible 
to those obtained with t y = 12; we only report results from the latter case. In order 
to extract the matrix elements and errors, we have performed a weighted average of 
R3 and X over t x . The intervals in t x were chosen on the basis of a study of the time- 
dependence of the effective energy (mass) obtained from G5 for a meson at rest. In 
order to have a good isolation of the lightest pseudoscalar state, and avoid exceedingly 
large statistical errors, the optimal time interval was found to be 10 < t x < 22. This 
was the case for both the Wilson and SW-Clover actions and for all the values of 
the hopping parameter. We estimate that the contamination from the excited states 
induces a systematic error of some per cent on the relevant quantities for all values 
of the quark masses considered in this study. Since as meson sources we used local 
operators, whose couplings to the physical states are independent of the momentum, the 
same interval in t x has been used for all meson momenta. Due to periodic boundary 
conditions, we also used the interval 42 < t x < 54. For R3, this corresponds to a 
different matrix element. In the following, X and R3 will always denote the weighted 
averages on t x . As mentioned above, computer memory limitations have obliged us 
to use the so-called "thinning" approximation. Thinning has been implemented here 
exactly as in p3] (which may be consulted for details). Although we have explicitly 
tested that thinning does not influence G5 , we have no way of testing its effect on more 
complicated correlations, such as Gq. We strongly suspect that this approximation 
results in an appreciable worsening of the statistical quality of our data. The rather 
large errors are also due to the relatively small size of our lattice. 
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4 Chiral behaviour of the Bk parameter 



We now examine the chiral behaviour of the matrix element, which we parametrize as 

(K°\d AS=2 \K°) = a + (3m 2 K + 7 (p • q) + 

+5m% + em 2 K (p ■ q) + ((p ■ qf + . . . (11) 

where all the quantities are expressed in lattice units and the ellipses stand for higher- 
order terms in (p ■ q) and rr? K . The matrix element on the l.h.s. of eq. ( |TT] ) is obtained 
by multiplying R3 of eq. ( pToj ) by Z5, which is extracted in the standard way from the 
asymptotic behaviour of the pseudoscalar two-point correlation G5 of eq. (||). The 
parameters a and (3 are lattice artefacts that should vanish in the continuum limit. 
As we have already stressed in section ^ their origin may be attributed to two effects: 
(i) discretization errors of the lattice matrix elements arising from the subtraction of 
operators with the wrong "naive" chirality; (ii) higher-order perturbative corrections 
to the renormalization constants Z + , Z$p, Zva, Z$pt of eq. @). Similar artefacts also 
contaminate the other parameters 7, 5, e, We investigate the differences in the values 
of a and /3, obtained by fitting the dependence of the matrix element on the kaon 



masses and momenta. Since in the continuum a and /3 vanish [ 30 1 , we consider a 
reduction of their values as a measure of the correctness of the chiral behaviour and of 
the accuracy in the determination of the matrix element. 

There are several possibilities for extracting the parameters a, (5, . . . , £. We have 
verified that they all give compatible results. In agreement with the observation of 
ref. JjJ, we found that, in order to minimize higher-order terms in m 2 K , the most 
convenient procedure is to substitute the fitting variables as follows 



>»i - A" - ^^§2 ™ic (12) 

O Zj^Zj 

2 



(p.,) => Y -^-\J^fp. q ) (13) 

and fit directly the ratio R3 with respect to X and Y. Up to quadratic terms, with 
a, /?,..., C appropriately redefined, eq. (|Tl]) becomes 

#3 = — (K \6 AS=2 \K°) = a + 0X + -yY + 5X 2 + eXY + QY 2 . (14) 

Another advantage is that R3 (and X) are obtained directly from ratios of Green func- 
tions which have correlated fluctuations. This has to be contrasted with {K°\0 AS=2 \K°) 
which is obtained by multiplying i?3 by Z5, computed from a fit of G*,(t x \Q). The re- 
duction of the size of the higher-order terms in m 2 K , obtained by using X and Y, is 
reflected in the stability of the results for a, (3 and 7. The values obtained using eq. 
([l4]) differ by only about 20% from those obtained from the linear fit 

R 3 = J-(K \6 AS=2 \K ) = a + (3X + ^Y. (15) 
Z5 
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Using eq. (|i~l|) instead, the results, although compatible within the errors, would differ 



by about 50%. Even though eq. (14) more than halves the discrepancy between linear 
and quadratic fits, the remaining systematic effect due to higher-order corrections 
requires further investigation. Below, we only present results from eqs. (^)- The 
S-parameter is computed from the relation 

B K ( j i)= 1 /Z\. (16) 

In the remainder of this section, we discuss possible systematic effects on the chiral 
behaviour of the matrix element of AS=2 . 



4.1 Chiral behaviour in perturbation theory 

One possible source of systematic error is that the one-loop perturbative determination 
of the renormalization constants, at current values of /?, is not accurate enough to 
correct the chiral behaviour. Thus, it is essential to improve the accuracy of the 
mixing renormalization constants which are calculated in PT. For this reason, the by 
now standard procedure of BPT [^lj has been adopted for both actions; we have tried 
two different boosted schemes, BS1 and BS2, by using the following definitions of a s 
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BS1: = 7^gy« s ~ 1.68 a s (W); 1.68 a s (SW) (17) 

BS2: of) = {%K c fa s ~ 2.49 a s (W); 1.84 a s (SW). (18) 

As can be seen, the effective coupling depends quite significantly on the prescription 
in the Wilson case, while the variation in the Clover case is small. In the following, 
we will also give results with ai 3) = 3.1 a, (BS3) pT[ 5 in the Wilson case, in order 



to compare with || and will call SPT the case in which we use the bare lattice a s 
as expansion parameter. In fig. 0, we show all our results for R3, plotted against 
Y = X(p • q))/m 2 K , where X = {Sf^m 2 K ) / {?,Z^Z\). Y > corresponds to the matrix 
element (K°\6 AS=2 \K ), whereas Y < corresponds to (K°^°|6 A5=2 |0) (with the 
usual caveat of final state interaction [29] and finite size effects) []. Notice that the 
operator is renormalized in the MS DRED scheme in the Wilson case, and in the RI 
scheme in the Clover case. 

In the following, we will only analyse data related to (K°\6 AS - 2 \K ). In fig. g, 
the effective coupling of the BS1 (which is the same for the two quark actions) was 
used to define the renormalized operator. In table we give the values of a, (5 and 7 
obtained from a fit of our data to eq. (|i~5|). Figure ^ and the values of a, /3 and 7 in 
table allow several observations: 



1 Lattice artefacts are clearly visible in the points corresponding to the highest momenta on the leftmost 
side of the figure. 
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Figure 2: Dependence of (0\d AS=2 \K°K°) (Y < 0) and (K \O AS=2 \K°) (Y > 0) on mass 
and momentum obtained with (a) the Wilson action (200 configurations) and (b) the SW- 
Clover action (460 configurations). The renormalized operators have been obtained by using 
oSp as the boosted coupling. 



Renormalization 


a 


P 


7 


SPT-W 
BSl-W 
BS2-W 
BS3-W 


-0.100(7) 
-0.068(5) 
-0.032(2) 
-0.0082(6) 


0.21(8) 
0.15(5) 
0.08(3) 
0.020(7) 


0.56(7) 
0.40(5) 
0.21(3) 
0.057(7) 


SPT-SW 
BS1-SW 
BS2-SW 


-0.067(12) 
-0.054(12) 
-0.052(12) 


0.17(15) 
0.17(15) 
0.16(15) 


0.62(11) 
0.62(11) 
0.62(11) 


NPM-SW n 2 a 2 = 0.96 
NPM-SW /i 2 a 2 = 2.47 


0.017(13) 
0.022(14) 


0.21(17) 
0.23(17) 


0.70(12) 
0.72(12) 



Table 2: SW-Clover (SW) and Wilson (W) results with different choices of mixing renor- 
malization constants. For comparison, in the Clover case the non-perturbative results at two 
values of the renormalization scale are also given. 
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1. The use of the SW-Clover action by itself does not have a significant effect on the 
chiral behaviour of the matrix element. 

2. For all the choices of boosted coupling, the fitted values of a in table |2| are different 
from zero. 

3. A comparison of boosted and unboosted results shows that, in the Clover case, a, 
(3 and 7 are almost unaffected by the boosting procedure. This is to be contrasted 
by the results in the Wilson case, which are rather unstable and strongly depen- 
dent on the precise choice of the effective coupling constant^ This is particularly 
relevant for 7, from which we estimate the Bk parameter. It implies that the 
estimate of Bk in the Wilson case is subject to a large systematic uncertainty, 
due to its sensitivity to the choice of the expansion parameter, i.e. to higher-order 
effects. 

4. In the Wilson case, the instability is due to the large value of F + (see table [l]): 
a, (3 and 7 diminish rapidly with increasing a s , roughly proportionally to Z + . 

Inspired by mean field theory we have used eq. (|j), with Z + as an overall 
factor, also in the perturbative case. Alternatively, we can use the formula 

A5=2 (/i) = Z + ( f ia,a s )d AS=2 (a) + Z SP 6sP (a) + Z VA V A (a) + Z SPT d S PT (a) , 

(19) 

which is equivalent at one loop, but may lead to significant differences for the 
matrix elements at large values of a s . Indeed, in the BS3 the two possibilities 
correspond to values of 7 which differ by a factor of 2.5. 



4.2 Chiral behaviour with the non-perturbative renor- 
malization 

We now examine the case in which the renormalization constants are determined non- 
perturbatively. This will only be done in the Clover case, where the mixing coefficients 
have already been computed. The data for the operator renormalized at fj, 2 a 2 = 0.96 
are shown in fig. ||[ In ref. |23j, the non-perturbative method of ref. [24] was imple- 



mented for the calculation of Z + , Z$p, Zva and Z$pt- For completeness, we reproduce 
the results of ref. |2^] in table [3|. With the NPM, the renormalization conditions for 

qAS=2 

can be imposed at any arbitrary scale fJ. 2 a 2 , subject to the conditions that lat- 
tice artefacts are small and that [i ^> Aqcd , so that continuum perturbation theory can 
be applied. In p3j, a window of values around fi 2 a 2 ~ 1 was established, for which the 
non-perturbative estimates of the Z's are fairly stable, see table HI the results in table ||] 
show that, with the non-perturbative renormalization, a and are always compatible 
with zero and 7 is very stable (for /i 2 a 2 > 0.66) in all the range of scales considered in 
our study. 



2 Notice that in the literature boosting factors varying between 1.6 and 3.1 have been used for this 
at /3 = 6.0. 
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u 2 a 2 


z+ 


L %sp 


Zva 


ZsPT 


0.46 


0.91 ±0.05 


0.08 ±0.14 


0.34 ±0.03 


0.34 ± 0.07 


0.66 


0.84 ± 0.05 


0.07 ±0.09 


0.33 ± 0.03 


0.34 ± 0.06 


0.81 


0.83 ± 0.04 


0.11 ±0.07 


0.31 ±0.03 


0.28 ±0.05 


0.96 


0.84 ± 0.03 


0.14 ±0.07 


0.30 ± 0.02 


0.24 ± 0.04 


1.27 


0.80 ± 0.04 


0.17 ±0.05 


0.29 ± 0.02 


0.21 ±0.03 


1.54 


0.82 ± 0.02 


0.19 ±0.04 


0.27 ±0.02 


0.16 ±0.03 


1.89 


0.83 ± 0.03 


0.22 ±0.05 


0.30 ±0.02 


0.18 ±0.03 


2.47 


0.85 ± 0.02 


0.22 ± 0.06 


0.33 ±0.02 


0.23 ± 0.03 


SPT 


0.91 


0.12 


0.12 


0.12 


BPT 


0.84 


0.21 


0.21 


0.21 



Table 3: Values of the Z's for several renormalization scales y?a 2 . We also give the re- 
sults obtained at n 2 a 2 = 1, by using "standard" perturbation theory (SPT) and "boosted" 
perturbation theory (BPT) with the BS1 effective coupling «W. 
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Figure 3: Dependence of {K°\0 AS=2 \K°) on mass and momentum, with the non-perturbative 
renormalization constants computed at ^i 2 a 2 = 0.96. The data can be compared with the 
corresponding perturbative ones in fig. 2. 
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The use of the values of the Z's in table [3| is subject at present to some limitations 
that we briefly want to discuss^. 

1. The calculation of ref. p3| was performed on 36 configurations at (3 = 6.0 on a 
16 3 x 32 lattice, using the SW-Clover action and "improved-improved" operators. 
The renormalization constants were determined for a single quark mass only, 
corresponding to K = 0.1425. We used these constants at all values of the quark 
masses. Even though for light quarks we expect very little dependence of the Z's 
on the quark mass [18, 28], one has to bear in mind that this approximation may 
introduce a small systematic effect in our results. 

2. The optimal strategy to obtain an accurate determination of the matrix element 
of the renormalized operator would be to compute the renormalization constants 
and the operator matrix elements on the same set of configurations and for the 
same values of the quark masses. This is not possible at the present stage, since 
the Z's and the matrix elements were computed on different sets of configurations, 
quark masses and lattice volumes. On the other hand, estimating the uncertainty 
of the final result by varying the Z's and the matrix elements independently within 
their errors is questionable because the values of the Z's, obtained by solving a 
linear system of equations on the same configurations, are strongly correlated. 
For this reason, we used the central values of the Z's from p3] and ignored 
their statistical error in the present analysis. To monitor the stability of the 
results, we varied over a wide range the renormalization scale at which the Z's 
were computed. Correspondingly, the Z's change by an amount comparable to 
their error, or even larger. A large set of results, obtained from a fit of the data 
to eq. fll5|), corresponding to 0.19 < fj?a 2 < 2.47 is given in table |3]. Notice 
that for scales below 0.81, the values of the Z's become unstable and the use of 
perturbation theory (for matching) questionable; for /j, 2 a 2 > 2.47 lattice artefacts 
become clearly visible |23| , 24 1. 

3. The stability of the results for the matrix elements of the renormalized operator 
appears surprising, since the renormalization constants vary appreciably on the 
same range of fi 2 a 2 [23]. The explanation for this is the following. The most un- 
stable of the renormalization constants, according to ref. [23], is Z$p. The Osp 
matrix element carries, however, the least weight in the final result. On the other 
hand, the most stable of the Z's is Zva, which gives the largest contribution. The 
case of the operator Ospt is intermediate between the two cases discussed before. 
We show this explicitly in table [5], where we present separately, for different quark 
masses, the contribution of the different operators to the final result. For compar- 
ison, we also give the different contributions in the perturbative case in the BS1. 
In the table, the contribution to R% of the bare operator O as=2 , corresponding to 
{K°\d AS=2 \K°)/Z 5 , will be denoted by R AS=2 ; the contribution of the operator 
Osp, corresponding to (K°\Osp\K } /Z5, etc., will be denoted by Rsp, etc. We 



These limitations are not intrinsic to the method and will be easily eliminated in future calculations. 
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/i 2 a 2 


a 


P 


7 


0.19 


0.03(2) 


0.2(2) 


0.99(17) 


0.46 


0.022(16) 


0.23(19) 


0.78(13) 


0.66 


0.013(15) 


0.20(17) 


0.70(12) 


0.81 


0.012(13) 


0.21(17) 


0.69(12) 


0.96 


0.017(13) 


0.21(17) 


0.70(12) 


1.27 


0.015(13) 


0.21(16) 


0.66(11) 


1.54 


0.018(13) 


0.22(16) 


0.67(12) 


1.89 


0.023(13) 


0.22(16) 


0.69(12) 


2.47 


0.022(14) 


0.23(17) 


0.72(12) 



Table 4: Values of the parameters of the fit of (K°\d AS - 2 \K°) to eq. (15). 



Method 


K 


R AS 2 


Rsp 


Rva 


Rspt 


Wilson 
BS1 


0.1550 
0.1540 
0.1530 


-0.004(5) 
0.027(4) 
0.058(4) 


-0.0154(4) 
-0.0161(3) 
-0.0165(2) 


0.0575(13) 
0.0620(10) 
0.0654(8) 


-0.0313(7) 
-0.0331(5) 
-0.0345(4) 


Clover 
BS1 


0.1440 
0.1432 
0.1425 


-0.024(5) 
-0.002(5) 
0.020(4) 


-0.0491(8) 
-0.0515(7) 
-0.0529(6) 


0.178(3) 
0.190(3) 
0.199(3) 


-0.1026(15) 
-0.1083(13) 
-0.1123(12) 


NPM-Clover 
H 2 a 2 = 0.96 


0.1440 
0.1432 
0.1425 


-0.024(5) 
-0.002(5) 
0.020(4) 


-0.0332(6) 
-0.0348(5) 
-0.0358(4) 


0.261(5) 
0.279(4) 
0.292(4) 


-0.1202(17) 
-0.1315(14) 
-0.1315(14) 


NPM-Clover 
/iV = 2.47 


0.1440 
0.1432 
0.1425 


-0.025(5) 
-0.002(5) 
0.021(4) 


-0.0540(9) 
-0.0567(8) 
-0.0582(7) 


0.284(5) 
0.305(5) 
0.319(4) 


-0.1161(17) 
-0.1225(15) 
-0.1270(14) 



Table 5: Separate contributions to R3 from the bare operator R AS=2 and the off-diagonal ma- 
trix elements, Rsp, Rva and Rspt, for different quark masses. We give the non-perturbative 
results at two different renormalization scales, and the perturbative result in the BS1. All 
the matrix elements correspond to the case where the two mesons are at rest. 
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see that the largest subtraction comes from the operator Ova, which is multiplied 
by the constant ZyA, which is quite well determined. This explains the stability 
of the fitting parameters a, (5 and 7 with varying renormalization scale. 

4. The improvement in the chiral behaviour of the matrix element arises from the 
difference in Zsp : ZyA ■ Zspt between the non-perturbative and perturbative 
case. Since ZyA is much larger than Z$p and Zspt in the non-perturbative case, 
most of the beneficial effect on the chiral behaviour is due to the enhanced contri- 
bution of Ova- Boosted perturbation theory, on the other hand, can only change 
all the coefficients by the same amount, leaving Z$p ■ ZyA ■ Zspt unaltered, 
so that the correction is much less effective. A very large boost of the coupling 
would be necessary in order to obtain the same effect as in the non-perturbative 
case. 



5 Physics results 



In order to compare the results for the renormalized ^-parameter obtained from 7, cf. 
eq. (|l~6| ) , with other theoretical predictions of the same quantity, it is useful to refer to 
a "standard" definition, which we will take to be the renormalization group invariant 
-B-parameter Bk- In the standard perturbative approach, by a suitable choice of Z+, 
the i?-parameter is defined in the MS scheme at a scale /i ~ 1/a. Let us call this 
i?-parameter B MS (fj 1 ). At the next-to-leading order (NLO), Bk is written in terms of 
B m {ji) as 



Br = ot s (ji)' 



7 (0) /2A) 



a s (n) / 7 (1) A) - 7 {0) A 



47T 



2/3 2 



(20) 



where /5o,i and 7( 0,1 ) are the leading and next-to-leading coefficients of the /3-function 
and anomalous dimension. /?o,i and 7^ are universal while depends on the renor- 
malization scheme. The explicit expressions of all the quantities appearing in eq. (p0|), 
^(oM 1 ) anc [ /3 0j i, can be found for example in refs. [32, |33f| . Notice that the renormal- 
ization group invariant i?-parameter defined in this way is also regularization-scheme 
independent, up to next-to-next-to-leading order terms. 

We compute Bk only for the Clover case using the NPM, since this is our best result. 
By using the non-perturbative estimate of Za = 1.06|1£], we obtain the S-parameter 
corresponding to the operator renormalized in the Regularization Independent (RI) 
scheme |2|, The matching to the MS scheme can be done in continuum perturba- 
tion theory (this is one of the main advantages of the non-perturbative approach), by 
computing the operator matrix element in the same gauge and on the same external 
quark states as those used for the non-perturbative calculation of the lattice renor- 
malization constants [23, p4fl. At the next-to-leading order, the NDR-RI matching 
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coefficient for O as 2 is given by B 



B MS (») = (l + ^Arf^j BjtO*) (21) 

where 

Ar MS = _ 14 / 3 + g log ( 2 ) _ ( 2 2) 

In numerical estimates, we have used a s (fj,) = 0.3, at the renormalization scale y?a? = 
0.96, corresponding to fj, ~ 2 GeV. Given the still large error on we have not 
included the uncertainty due to the choice of a s (fj,) in the result given below. At fixed 
a s (/j,), there is still a variation of ±0.03 when varying rif, the number of active flavours, 
between and 4 in eq. (pC| ). With the above parameters, we get 

Br = 0.86 ± 0.15 . (23) 

Within large errors, our results agree with other lattice calculations and with other 
results obtained with a different theoretical approach, namely the 1/N C expansion 



p% pq| . It may be noted that, in spite of our increased statistics, the final errors 
for our results are comparable (or even larger) to those of previous lattice studies. 
We attribute this to our thinning approximation, which we plan to remove in future 
calculations. 



6 Conclusions 

We have presented a parallel lattice calculation of the kaon B parameter Bk with the 
Wilson and the SW-Clover actions, at (3 = 6.0. In the Wilson case, using perturbation 
theory to determine the renormalized four-fermion operator, we found that the value 
of Bk strongly depends on the precise choice of the expansion parameter a s and 
on higher-order corrections. Consequently the results are subject to a rather large 
uncertainty. In the Clover case, the results were more stable with respect to a change 
of a s . The most appreciable improvement in the chiral behaviour of the operator is, 
however, found with the operator renormalized non-perturbatively, by using the mixing 
coefficients computed in ref. [23]. 



The accuracy of our results was somehow limited by two reasons: the precision 
in the calculation of the mixing coefficients and the memory of the APE 6-Gflops 
machine, which forced us to work on a relatively small lattice volume and to use the 
"thinning" approximation. These limitations are not intrinsic to the non-perturbative 
method and will be easily eliminated in future calculations. In particular, we expect 
significant improvements by computing renormalization constants and operator matrix 
elements on the same set of configurations and for the same values of the quark masses 
(and, of course, by eliminating the "thinning" and by working on a larger volume). 
We also plan to compute non-perturbatively the mixing coefficients in the Wilson case, 
although, in this case, discretization errors could still play an important role. In the 
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non-perturbative case, with more accurate results, it will be possible to check whether 
a and f3 are so small as to give a negligible contribution to the matrix element for 
quark masses corresponding to the kaon. Then, similarly to the case of staggered 
fermions, it will be possible to compute the kaon matrix element directly, without 
the fitting procedures used so far. This will allow us then to reduce the error in the 
determination of the quenched Bk to about 5%, as in the staggered case. 

The results of our study are encouraging, and motivate us to extend the calcula- 
tion of matrix elements of operators renormalized non-perturbatively to the operators 
relevant to Al =1/2 transitions and to the penguin and electro-penguin operators 
which control CP-violation in kaon systems. Our strategy is to achieve an accurate 
determination of the physical weak amplitudes, by combining the improvement of the 
action a la Symanzik, which reduces 0(a) effects, with the non-perturbative method 
ofref. @. 
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